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We perform a numerical analysis of Higgs-to-Higgs decays within a Type II 2-Higgs Doublet Model 
(2HDM), highlighting several channels that cannot occur in its Supersymmetric version, thereby 
allowing one to possibly distinguish between these two scenarios. Our results are compliant with 
all available experimental bounds from both direct and indirect Higgs searches and with theoretical 
constraints from vacuum stability and perturbative unitarity. 



I. INTRODUCTION 

Whilst Supersymmetry (SUSY) [l| provides an attractive theoretical scenario for physics beyond the 
Standard Model (SM) (in its ability to remedy the hierarchy problem, to provide a natural dark matter 
candidate, to enable high scale gauge coupling unification, etc.), there is to date no evidence for it. In 
its miminal (in terms of particle content and gauge structure) incarnation, the Minimal Supersymmetric 
Standard Model (MSSM) 0, SUSY is highly predictive though, e.g., in the Higgs sector. Of the initial 
eight degrees of freedom pertaining to the two complex Higgs doublets responsible for Electro- Weak 
Symmetry Breaking (EWSB) in the MSSM, after the latter has taken place, three are consumed to give 
mass to the SU{2) g> U{\) weak gauge bosons, and Z, so that five physical Higgs states survive: 
the CP-neutral ones h and H (Mh > Mh), the CP-odd one A and the charged states H^. In effect, 
regarding the Higgs sectors the MSSM is nothing but a 2-Higgs Doublet Model (2HDM) of Type II (in 
the nomenclature of Ref. [3J]), whereby SUSY enforces relations amongst Higgs masses and couplings, on 
the one hand, and weak gauge boson masses and interaction parameters, on the other hand, in such a 
way that, of the original 7 independent parameters defining a CP-conserving 2HDM (see below), only two 
survive as such in the MSSM [J|. These can be taken to be tan/3, the ratio of the Vacuum Expectation 
Values (VEVs) of the two Higgs doublet fields, and one of the extra Higgs boson masses. At tree-level, 
these are the only inputs needed to compute Higgs masses and couplings to ordinary matter (quarks, 
leptons and gauge bosons) as well as Higgs self-couplings. Hence, if SUSY had chosen all the sparticle 
masses to be much larger than the SM objects and the Higgs bosons, so that they are not accessible 
at the upcoming Large Hadron Collider (LHC), one intriguing question to ask would be whether it is 
possible to distinguish between the MSSM and a Type II 2HDM on the sole basis of Higgs interactions 
with SM matter and/or self- interactions. Or, conversely, whether it would be possible to dismiss the 
assumption of such (decoupled) heavy sparticles, hence of minimal SUSY altogether (aka the MSSM), 
from the observation of particular signals in the Higgs sector alone. 

It is the purpose of this work to prove that this is the case, exploiting the fact that SUSY prevents 
some Higgs-to-Higgs decays, that remain instead possible in a generic Type II 2HDM. This is primarily 
connected to the fact that a generic pattern of Higgs masses, as dictated by SUSY, in the MSSM is the 
one in which the h is rather light (in fact, below 130 GeV or so, for heavy sparticles [5() whilst the others 
(H, A and H ± ) are quite heavy and degenerate in mass, the more so the larger tan/3 (33|. In fact, even 
in the presence of off-shellness effects [g, 0] in all decay products in Higgs decay chains in the MSSM, 
essentially only the H — > hh and A — > Zh channels are possible in this scenario [34]. In contrast, in a 
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Type II 2HDM without SUSY, many other decays are possible, e.g.: 

H — > AA, H — > H + H~ , H^W ± HT : H ^ ZA, 
A->W ± H^ 1 A->ZH, 
-> W ± h, H ± -> W ± H -> W^A (1) 

While the existence of such different decay patterns in the two models has been known for some time [|| , 
our ultimate intention here is to prove that the 2HDM Type II is still phenomenologically viable in the 
light of all available experimental constraints. In practise then, whilst, of course, these decays could not all 
occur at the same time, depending on the Type II 2HDM parameters, a subset of them would be possible, 
therefore providing a means of distinguishing between the two scenarios discussed, further considering 
that - under the above assumption of a heavy SUSY spectrum - the dominant production channels of 
both neutral and charged Higgs bosons in both scenarios are the same and proceed via interactions with 
SM particles 

The plan of the paper is as follows. In the next section we describe the Type II 2HDM that we will be using 
and we fix our conventions. The following section discusses the experimental and theoretical bounds. We 
then illustrate our results quantitatively in section IV. The final section outlines the conclusions. 



II. THE TWO-HIGGS DOUBLET MODEL 



We start with a brief review of the 2HDM used in this work. The potential chosen is the most general, 
renormalisable and invariant under SU(2) <8> U(l) that one can build with two complex Higgs doublets 
with a softly broken Zi symmetry. It can be written as 



U 2 hdm = A*?|$i| +l4\®2 



(nl<S>\§2 +h.c.) 
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where $i and $2 are the two Higgs doublets with hypercharge +1/2 and /x§ is the Z2 soft breaking term. 
For simplicity we can take /i§ and A5 to be real. The doublet fields are parameterised as 



+hi + izi) 



(i = 1,2), 



(3) 



(246 GeV) 2 . Assuming CP-conservation, this potential 



where the VEVs v\ and v-i satisfy v\ + v\ = 1 
has 8 independent parameters. However, because v is fixed by the boson mass, only 7 independent 
parameters remain to be chosen, which we take to be Mh, Mjj, Ma, M h ±, tan/3, a (the mixing angle 
between the two CP-even neutral Higgs states) and M 2 — /i 2 /(sin (3 cos (3) (which is a measure of how the 
discrete symmetry is broken). The definition of a and (3 and the relation among physical scalar masses 
and coupling constants are shown in Ref. [l(| for definiteness. 

In a general 2HDM, the Yukawa Lagrangian can be built in four different and independent ways so that 
it is free from Flavour Changing Neutral Currents (FCNCs). We define as Type II the model where <f>2 
couples to up-type quarks and <t>\ couples to down-type quarks and leptons [35| . (We present the Yukawa 
couplings for a Type II 2HDM in the Appendix.) 



III. EXPERIMENTAL AND THEORETICAL BOUNDS 



The results from all LEP collaborations on topological searches with two Higgs bosons or one Higgs and 
one gauge boson were presented in Ref. [l2j]. We will use the experimental limits on the cross sections 
for e + e~ — * H\Z and e + e~ — » H1H2, where Hi can be any CP-even Higgs boson and H2 can be either 
a CP-even or a CP-odd Higgs boson. As we are concerned here with a Type II 2HDM, there is a bound 
particularly relevant to our analysis which is the one obtained from the relation involving the following 
cross sections (a) and Branching Ratios (BRs): 

aye e — * n 1 A) , „2HDM u~u\ ■ 2/ a\ r,T3 ( rr2HDM j7\ , ,\ 

a( e + e - _> Hf M Z) BR ^ Hl hb > = 8m ^ a ~ BR ( H i ~* bb > ■ ( 4 ) 
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In a Type II 2HDM h — > 66 is the dominant decay for most of the parameter space for Mh below the SM 
limit. The two subleading competing decays are h — ► cc and h -> t+t~. In a Type II 2HDM, the ratio 
T(h — > bb)/T(h — > t + t~) is the SM one. In contrast, one has 



r(/t -> 6b) 

r(/i — > cc) 



'6 



tan 2 a tan 2 (3 (5) 



in the limit mh >> m q . If we then take the case a rj /3 and because we always consider tan/3 > 1 the 
decay /i — > bb becomes even more dominant in a Type II 2HDM, with respect to the SM case. Conversely, 
if a w [3 + 7r/2, then we recover the SM ratio. Either way, h — ► 66 is the dominant decay by an amount 
which is at least the corresponding SM ratio of BRs with respect to the other fermionic decays. Because 
its dependence on the other parameters of the model is very mild, this essentially provides a limit in the 
(3 — a versus Mh plane. In particular, it is straightforward to check that when sin(/3 — a) rs 0.1 there is 
essentially no bound on the lightest Higgs boson mass, while for sin(/3 — a) « 0.2 the limit immediately 
jumps to Mh > 75.6 GeV. 

In this study the masses of the charged Higgs boson and the CP-odd one will always be above 200 GeV 
and therefore not constrained at all by the LEP bounds from direct searches. In some cases though, the 
mass of the heavy CP-even Higgs boson will be allowed to be below 200 GeV. In those instances we have 
confirmed that the LEP bounds do not apply to the cases presented in this work. Concerning if ± states, 
apart from a model independent LEP bound of Aljj± > My/± , DO and CDF have model dependent limits 
on the charged Higgs mass (see Ref. [l3j) from top decays, but again these are below the range of H ± 
masses discussed in this work. No other experimental bounds exist from direct searches for the set of 
parameters that we will present. 

Other than limits from direct searches for Higgs bosons, there are indirect constraints from precision 
observables, from both LEP and SLC. New contributions to the p parameter stemming from Higgs states 
[13 ] have to comply with the current limits from precision measurements [l5j]: \5p\ < 10 -3 . There are 
limiting cases though, related to an underlying custodial symmetry, where the extra contributions to dp 
vanish. In this study we will consider two such particular cases: (A) the one where M H± = M A and 
(B) the one where M H± = M H with sin(/3 — a) = 1. These parameter choices correspond to the case 
in which the custodial symmetry (SU(2) L <8 SU(2) R — ► SU{2) V ) is preserved in the Higgs potential, so 
that the latter can be written in terms of Tr(MjM/) with M l — (ir2<&|, where the Mi's (i — 1,2) 
are translated as Mj — > M[ — g[Mig R with g L R £ SU{2) L . R (A) or Tr(M 2 iM 2 t 1 ) and det(M 2 i), with 
M 2 i = (iT 2 $27*i) ; where M 21 is translated as M 2 i -> M' 2X = g\M 2 \g R with g LR E SU{2) LtR (B), 
respectively [T3 |. In Fig. [U we show the allowed region under the p parameter constraint in several 
scenarios which are relevant to our later discussions. Furthermore, it has recently been shown in Ref. [l7| 
that, for a Type II 2HDM, data on B — > X s j imposes a lower limit of M H ± > 250 GeV, which is 
essentially tan/3 independent. Other experimental constraints on a Type II 2HDM come from the results 
on (g — 2) u (the muon anomalous magnetic moment) jlal . Rb (the 6-jet fraction in e + e~ — > Z — > 
jets) [l^, H3|, the decay B + — > t + v [2l| . B q B q mixing [I3| and the r leptonic decay In general, 

bounds from these observables can be important for relatively small values of M H ± and large tan /3 > 10 
[l8l [T^ . [20l . |2ll [22I [23l . [2^ | . Values of tan (3 smaller than « 1 are disallowed both by the constraints coming 
from Z — > 66 and from B q B q mixing. 

Concerning theoretical constraints we will take all masses Mh, Mh, Ma and M#± to be below 700 GeV. 
This is a consequence of tree-level unitarity bounds (2f| [2(| in the limit M = 0. Furthermore, the most 
general set of conditions for the Higgs potential to be bounded from below are [13] 

Ai>0, A 2 >0, y/\! A 2 + A 3 + min (0, A 4 - I A 5 1) > 0. (6) 

Recently, it was in fact proven that these are necessary and sufficient conditions to assure vacuum 
stability of the potential at tree level [28|. Vacuum stability against charge breaking is also built into a 
Type II 2HDM model, as a non-charge breaking minimum, when it exists, is always the global one in any 
2HDM [2^|. Finally, according to [3(J, perturbativity for the top and bottom Yukawa couplings forces 
tan/3 to lie in the range 0.3 < tan/3 < 100, though it turns out that, from enforcing perturbativity also 
on the Aj's, moderate values of tan/3 (tan/3 ~ 0(1)) are preferred, especially for M ~ GeV. 

In the following, we will consider M H ± to be 250 GeV or larger and take tan/3 ~ 1-3, so that all the 
Type II 2HDM scenarios presented are free from these bounds. We will show that even with such small 
values of tan/3 a lot of the parameter space is already excluded. This does not mean that larger values 
of tan/3 are not allowed but rather that they are less likely to occur. Choosing a large tan/3 forces a 
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FIG. 1: Allowed regions under the p parameter data (at 2<r level). [Left figure] The area surrounded by solid 
curves corresponds to the allowed regions for the case of Mh = 120 GeV, M H ± = 250 GeV, sin(/3 — a) = 1. 
The one surrounded by dashed curves corresponds to the allowed regions for the case of Mh = 50 GeV, M H± — 
250 GeV, sin(/9 — <a)=0.1. [Right figure] The area surrounded by solid curves corresponds to the allowed regions for 
the case of Mh = 120 GeV, M H — 150 GeV, sin(/3 — a)=0.9. The one surrounded by dashed curves corresponds 
to the allowed regions for the case of M h = 50 GeV,M H = 150 GeV, sin(/3 — a)=0.1. 

very particular set of values for the remaining free parameters if one is to comply with all constraints. 
Therefore it may seem we are just scanning over small corners of the 2HDM parameter space. This is not 
the case. The 2HDM is already tightly constrained by experiment and it becomes severely constrained 
when one adds the theoretical bounds, especially those from perturbative unitarity. The EW precision 
data are also very restrictive. As shown in Fig. [T] (left), where we have fixed the mass of the charged 
Higgs at 250 GeV, a vast range of values of the heavy CP-even Higgs is possible. The allowed region 
for the pseudo-scalar masses depends on the values of (3 — a: when sin(/? — a) ~ 1 (full curve) all values 
of Ma are allowed provided Mh is close to M H ± = 250 GeV. As we move away from sin(/3 — a) = 1 
(dashed curve) the range of allowed values of the pseudo-scalar mass shrinks to a region around the value 
of M H ± — 250 GeV as wide as the precision measurements permit. On the right plot we see exactly the 
same trend but now in the M H ± vs Ma plane. In the plots shown in the following sections we choose the 
exact limits that cancel the p parameter contribution. We note however that we have explicitly checked 
that varying the values of the masses complying with these constraints does not produce a qualitative 
change in our analysis. In most cases not even a quantitative change is noticed. We end this section 
by underlining once more that we are not focusing on a small corner of the 2HDM parameter space. 
Considering all constraints a 2HDM Type II is subject to and the above discussion it is clear that we are 
spanning the entire parameter range allowed. When we focus on a definite limit, like M H ± = Ma, it is 
for illustrative purposes only. 



IV. DECAYS 

Let us start by saying that all the widths and BRs presented in this work are calculated at tree-level except 
for the decays to gg (plus 77 and Z7, not visible though in our plots) which are one-loop processes at the 
lowest order. However, we take into account the leading one-loop QCD corrections to Higgs to quark- 
antiquark decays by computing one-loop running masses for the (anti)quarks in the Yukawa couplings, 
evaluated at the decaying Higgs mass. 
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A. A decays 
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FIG. 2: A decays for M h = 120 GeV, M H = M H ± = 250 GeV and sin(/3 - a) = 1. On the left is the plot for 
tan/3 = 1 while on the right we set tan/3 = 3. We show the perturbative unitarity limits for M = GeV and 
M = 250 GeV on the left plot and for M = 200 GeV and M = 250 GeV on the right plot. For M = GeV 
and tan/3 = 3 all parameter space is excluded. We also show the vacuum stability limits for M — 250 GeV. For 
M = 200 GeV and below all parameter space is allowed in what concerns vacuum stability. 

For the CP-odd scalar and in the mass region chosen, sin(/3 — a) has to be very close to 1, in order for 
the model to be consistent with experimental data. Even a value of sin(/3 — a) = 0.9 is enough to violate 
precision measurements via the p parameter. Therefore the decay A — > Zh is not allowed. We take 
Mh = 120 GeV but the CP-odd Higgs boson profile does not depend on the light CP-even Higgs boson 
mass. We have chosen Mjj = M H ± — 250 GeV due to the B — > X s j bound and to the p constraint. 
Fig. [2] illustrates the decay patterns of the CP-odd Higgs state for two values of tan (3. Two comments 
are in order here. Firstly, notice that to choose larger values for Mjj and M H ± would only have the effect 
that the corresponding channels would open later. Secondly, the dependence on tan (3 is generally as 
follows: the larger tan/3 the more suppressed the decay into tt and consequently the CP-odd Higgs state 
decays more and more into other Higgs bosons as soon as the corresponding channels are kinematically 
allowed. All decay channels shown in these plots do not depend on M, as no Higgs self-couplings are 
involved in those processes. However, both the perturbative unitarity and the vacuum stability bounds 
depend on the value chosen for M. The excluded regions due to the the above constraints for the three 
values of M = 0, 200 and 250 GeV are shown in Fig. O Note that the smaller the parameter M is the 
smaller tan/3 has to be to avoid the perturbative unitarity limit. On the contrary, the smaller M is, the 
less constrained is the parameter space from the vacuum stability conditions (for M — GeV no bounds 
apply). 

In the MSSM, for tan j3 = 3, the pseudo-scalar decays mainly to fermion pairs, bb in the low mass region 
and tt when the channel becomes kinematically allowed in most of the studied scenarios [3l|. There 
is a situation where the MSSM and the general 2HDM are hard to distinguish. The branching ratio of 
A — > ZH can be very similar in both models when we compare the so-called intermediate-coupling regime 
of the MSSM (tan/3 w 3 and H/A masses below the the tt threshold) with a 2HDM with a large charged 
Higgs mass so that the the decay to W + H~ is forbidden. In Fig. [3] we show the total A width for the two 
situations presented in Fig. O When the A decays mainly to bb or gg the width is negligible but, when 
the ti, the other Higgs or Higgs plus gauge boson channels open, Ta grows rapidly reaching 100 GeV for 
Ma — 700 GeV. Note that the theoretical constraints shown in Fig. [5] are not shown again in this plot. 
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FIG. 3: Total width of the CP-odd Higgs state A for for the parameter values presented in Fig. [2] 
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FIG. 4: H decays for M A = M H ± = 250 GeV and tan/3 = 1. On the left is the plot for M h = 50 GeV and 
sin(/3 — a) = 0.1 while on the right we have Mh — 120 GeV and sin(/3 — a) = 1. We show the perturbative 
unitarity and the vacuum stability limits, for M = OGeV for the left plot and for M = GeV and M = 250 GeV 
for the right plot. 

B. H decays 

In this section we deal with the decays of the heavier CP-even Higgs boson. New contributions to the p 
parameter are avoided by setting Ma = M H ± = 250 GeV. We note once more that, as shown in Fig. [TJ 
this limit can be relaxed. As a consequence the dominant channel can be either the one with a charged 
Higgs or the one with a pseudo-scalar boson depending on the relation between the respective masses. 
Again, the constraint from B — > X s j is used. Then, we distinguish between two extreme situations. The 
first one, shown in Fig. [4] (left), is for sin(/3 — a) = 0.1. In this case, the H couplings to the gauge bosons 
are close to the SM ones. Due to the small value of sin(/3 — a) the mass bound on the light Higgs can easily 
be evaded and we choose the mass Mh = 50 GeV (though smaller masses are also allowed). This is also 
the limit where the H couplings to W ± H zfl and ZA are very small. In this left plot we can distinguish 
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FIG. 5: Total width of the CP-even heavy Higgs state H for the parameter values presented in Fig. [4] 



two interesting regions with new physics signatures. For small H masses, H — > hh, W + W~ , ZZ can be 
the dominant decays. For large H masses, the decays to a pair of charged Higgs bosons and to a pair of 
CP-odd Higgs states dominate as soon as they are kinematically allowed (though for a small Mh interval). 
The larger tan /3 is the more dominant these decays become. We also present the perturbative unitarity 
limit for M = GeV. We still present the decays that are above the perturbative unitarity limit for two 
reasons. Firstly, because as M grows the allowed Mjj region also grows (although for, for example, for 
M = 200 GeV, the H — > hh channel is negligible and in the low H mass region H — > W + W~ always 
dominate) and the decay to two charged Higgs and/or two pseudo-scalars are then allowed over a much 
larger Mh interval. Secondly, because, had we chosen a lower value for Ma and Mjj± , the corresponding 
decays would have opened for lower Mjj values, hence well within the region allowed by perturbative 
unitarity. 

The other extreme situation is shown on the right hand side of Fig. [3] and occurs for sin(/3 — a) = 1. 
In this limit, the H couplings to the gauge bosons are exactly zero. Decays to the two light Higgs states 
are also suppressed but they could still play a role if the soft breaking parameter is different from zero. 
This is the limit where H couplings to W ± H T and ZA are largest. Again we show perturbative unitarity 
and vacuum stability constraints for this case. Notice that the choice of tan j3 = 1 is heavily imposed by 
these bounds: for tan/3 = 2 and for the same set of parameters shown in the plot, the heavy CP-even 
Higgs mass is forced to be below w 350 GeV. Even if M is raised to 200 GeV the bound only grows 
to Ri 400 GeV. In Fig. [5] we show the total H width for the two situations presented in Fig. 2] When 
Mh = 50GeV, left plot in Fig. [U the heavy Higgs is allowed to decay to other Higgs and gauge bosons 
and that is why the H width is SM-like for the same mass. In the other scenario the heavy CP-even 
Higgs is not allowed to decay to gauge bosons. It decays to two gluons and fcrmion pairs in the low mass 
region. Therefore the H width is much smaller. As soon as channels with Higgs and gauge boson open 
both widths converge to similar values. 

In the MSSM, for tan/3 > 1, the heavy Higgs decays mainly to fermion pairs bb and then tt in the 
decoupling regime. Outside this regime there are two particular cases where distinguishing between both 
models will be hard. The first one is the anti-decoupling regime (tan/3 > 10 and Ma < M™ ax ) where, 
if kinematically allowed, the H — > hh can be the dominant decay channel. Hence, a 2HDM Higgs as 
presented in the left plot of Fig. |4] can be mistaken by such a MSSM heavy Higgs. The same is true 
for the above described intermediate-coupling regime where Br(H — ► hh) reaches 60 % for a significant 
heavy Higgs mass region. For a detailed discussion see |3lj |. 



C. H decays 



This section is dedicated to charged Higgs boson decays. Again, we choose M H ± = Ma to avoid the 
constraints from the p parameter. Once more we distinguish between two extreme cases regarding the 
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value of sin(/3 — a), which is the parameter that regulates the coupling to other Higgses and gauge 
bosons. When sin(/3 — a) is such that the LEP bounds can be avoided there are mainly two competing 
decays for the allowed charged Higgs boson mass region: H + — » tb and H + — > W + h. In Fig. [6] we show 
the charged Higgs BRs for sin(/3 — a) = 0.1 and Mh = 50 GeV for two values of tan/3. It is clear that 
the decay H + — » is always important and becomes dominant for large values of tan/3. 




FIG. 6: decays for M h = 50 GeV~, M H = 150 GeV, M A = M H ± and sin(/3 - a) = 0.1. On the left is the plot 
for tan f3 = 1 while on the right we set tan f3 — 3. Perturbative unitarity limits are shown. 



The other extreme case, sin(/3 — a) — 0.9, is plotted in Fig. [71 In this case the coupling to the heavier 
CP-even Higgs boson becomes dominant relative to the light Higgs case and the decay H + — > W + H 
is now the leading one for large values of tan/3. The p constraint could alternatively be enforced by 
M H ± w Mh because sin(/3 — a) w 1 as can be seen from the left plot in Fig. [1] In that scenario the 
HW + final state would be replaced by H + — > AW + which is independent of the value chosen for tan/3. 
Again, we take Mh = 150 GeV, though this value has no bearing on the final result except when we 




FIG. 7: H ± decays for M h = 120 GeV, M H = 150 GeV, M A = M H ± and sin(/? - a) = 0.9. On the left is the 
plot for tan /3 = 1 while on the right we set tan (3 — 3. Perturbative unitarity limits are shown. 

are in a region of large tan/3 and large sin(/3 — a). However, the large tan/3 region is excluded by the 
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FIG. 8: Total width of the charged Higgs state H for the parameter values presented in Figs. [6] and [7] 



perturbative unitarity constraints and therefore we will not consider this scenario. In Fig. [5] we show the 
total iJ ± width for the situations presented in Figs. [6] and [7l It is clear from the plot that the width does 
not depend so much on the parameters as happened for the previous cases. This is due to the fact that 
all channel types are already allowed starting from Mjj± = 250 GeV . 

For the mass regions considered in the plots, all MSSM scenarios predict an almost 100 % decay to tb. As 
stated in the introduction, the H ± — > W ± h channel although kinematically possible in the MSSM, only 
occurs with sizable rates in a tan region which is already excluded by experimental data. Therefore to 
distinguish a charged Higgs from 2HDM we should look for sizeable final states with one W boson and 
some other scalar. 



V. CONCLUSIONS 



We have highlighted that in a Type II 2HDM there exist Higgs-to-Higgs decays which are prevented from 
occurring in its SUSY version, the MSSM, owing to the fact that the latter imposes stringent relations 
amongst the masses of the H, A and 7? ± states, so that they are degenerate in mass over most of the 
parameter space. As these modes typically involve decaying Higgs states that are rather heavy, they 
could be primary means available at the LHC (and much less so at the Tevatron) to dispell the MSSM 
hypothesis that assumes that the sparticle states are very heavy and beyond the kincmatical reach of the 
collider. An analysis of Higgs pair production in the same spirit is now also in progress (32l |. 



APPENDIX A: YUKAWA COUPLINGS OF A TYPE II 2HDM 



In this Appendix we present the Feynman rules for the Type II 2HDM Yukawa couplings. 
Hereafter, the label u refers to up-type quarks and neutrinos whilst d to down-type quarks 
and leptons. Also notice that the Goldstone bosons couple just like in the SM, so we 
do not report their fcrmionic interactions here. Finally, we define 7^ = (1 — 7s)/2 and 
7i? = (1 + 7s)/2. Using notation already introduced (apart from Vij being the Cabibbo-Kobayashi- 
Maskawa matrix element in the quark sector and equating to 1 in the lepton case), one has: 



10 



Ti n h- ig cos a m ddh- ig sina m,, 

2M W sin/3 U,jU,f(.. 2Afw cos/3 

UidjH + : ^ 9 Mw Vij [tan /? TO dj 7^ + cot /3 m Mi 7 L ] d t UjH~: ^ 9 Mw Vfi [tan /? m rft 7 L + cot /3 m Mj 7^] 
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